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INTRODUCTION 

The  infiltration  capacity  of  a  soil  as  it  changes  with  continuing 
rainfall  is  a  most  important  estimate  in  computing  the  hydrologic  performance 
of  watersheds.   Many  of  the  designs  developed  in  watershed  engineering  are 
based  upon  man's  ability  to  modify  the  infiltration  capacities  of  the 
various  soils  within  the  watershed.   The  hydrograph  of  runoff  is  greatly 
influenced  by  infiltration.   A  dependable  system  is  needed  for  associating 
the  infiltration  curve  with  characteristics  of  the  plant  cover  and  the  soil 
profile.   The  concepts  presented  in  this  paper  offer  some  new  avenues  of 
approach. 

2/ 
As  reviewed  by  Philip  (7)  —  ,  the  classic  equations  of  Horton  and 

Kostiakov  require  knowledge  of  the  initial  rate  of  infiltration  and  compute 

deceleration  as  a  function  of  time  to  derive  the  curve  of  infiltration.   A 

modest  error  in  the  estimate  of  the  initial  rate  or  in  its  deceleration  can 

result  in  serious  miscalculations  of  infiltration  amounts,  especially  for 

larger  periods  of  time.   Philip  (7)  developed  an  equation  with  time  as  the 

dependent  variable  (i.e.,  the  rate  of  infiltration  a  function  of  cumulative 

infiltration)  but  abandoned  it,  because  of  its  unwieldiness,  in  favor  of  an 

equation  for  computing  the  infiltration  rate  as  a  function  of  time  with 

sorption  and  permeability  as  parameters.   Sorptivity  is  defined  by  Philip  as 

the  capacity  to  absorb  or  desorb  liquid  by  capillarity. 

This  paper  proposes  a  concept  for  estimating  a  potential  volume  of 
infiltration  from  characteristics  of  the  soil  and  subsequently  providing  a 
shape  to  the  curve  of  progress  toward  this  total.   The  rate  of  infiltration 
is  herein  considered  to  be  a  function  of  the  remaining  volume  of  potential 
storage  above  the  confining  horizon  and  the  permeability  of  the  confining 
horizon.   Time  is  therefore  a  dependent  variable,  which  is  a  function  of  the 
volume  increment  and  the  associated  rate.   The  unwieldiness  encountered  by 
Philip  is  not  dispelled  but  the  technique  is  demonstrated  by  simple  sequences 
of  computations. 

An  array  of  prominent  soils  in  the  United  States  prepared  by 
G.  W.  Musgrave  (6)  provides  a  means  of  evaluating  the  constant  rate  of 
infiltration  after  the  soil  is  thoroughly  wetted.   Since  Musgrave  derived 
his  estimates  of  infiltration  from  watershed  hydrographs  wherein  quick-return 

— '  Director,  U.  S.  Hydrograph  Laboratory,  Northeast  Branch,  Soil  and 
Water  Conservation  Research  Division,  Agricultural  Research  Service, 
U.  S.  Department  of  Agriculture,  Beltsville,  Maryland. 

2J   Numbers  in  parentheses  refer  to  Literature  Cited  at  the  end  of  this 
report. 


flow  is  indistinguishable  from  surface  runoff,  his  estimates  of  final  rate 
of  infiltration  approximate  the  permeability  of  the  least  permeable  horizon 
reached.   There  remains,  however,  the  problem  of  estimating  the  volume  and 
rate  of  infiltration  preceding  this  constancy.   An  inf iltrometer  survey  of 
two  small  watersheds  in  Illinois  in  1940  (8),  although  not  designed  for 
this  purpose,  provides  some  data  suitable  for  preliminary  exploration  of 
this  aspect.   The  symbols  used  in  this  discussion  are  given  as  footnotes 
to  table  1. 


POTENTIAL  VOLUME  OF  INFILTRATION 

The  volume  of  infiltration,  which  occurred  before  the  rate  of  infil- 
tration became  constant,  was  determined  with  an  inf iltrometer  for  various 
antecedent  moisture  contents  on  plots  selected  to  represent  the  diversity 
of  soils  and  vegetation  found  on  the  Edwardsville,  111.,  watersheds. 
Figure  1  portrays  a  "run"  typical  of  those  made  on  6  -  by  12-foot  rectan- 
gular plots  with  a  type  "F"  inf iltrometer  applying  about  1.78  inches  of 
water  per  hour.   Data  in  table  1  are  from  plots  on  which:   (a)   Recurrent 
runs  were  made,  (b)  one  type  of  vegetation  predominated,  and  (c)  no 
unusual  condition,  such  as  return  flow  at  some  down-slope  point,  was 
observed.   The  volumes  of  infiltration  prior  to  the  time  of  a  constant  rate 
(F,  @,  fc),  given  in  the  next  to  the  last  column  of  table  1,  were  derived 
as  illustrated  in  figure  1. 

Although  the  detail  characteristics  of  soils  on  these  plots  are  not 
essential  to  the  portrayal  of  the  concepts  presented  here,  descriptions  of 
typical  profiles  are  as  follows: 

Plot  No.  2         Alma  Silt  Loam         Slope         9.62  percent 

A-  0"  to  7".   A  light  brownish-yellow  friable  silt 
loam,  with  a  rather  heavy  grayish  cast  especially 
noticeable  when  dry.   Single  grain  structure. 
No.  A2  horizon  is  present,  indicating  that  all  the 
surface  soil  remaining  has  been  mixed  together  by 
cultivation. 

Bl-  7"  to  10".   A  light  brownish-yellow  heavy  silt 
loam  with  slight  grayish  coatings  and  some  gray 
mottling. 

B2-  10"  to  20".   A  light  brownish-  to  reddish- 
yellow  silty  clay  loam  rather  heavily  mottled  with 
light  gray.   Slightly  angular  aggregates. 

C-  20"  to  40".   A  mottled  light  gray  and  reddish- 
yellow  heavy  silt  to  silty  clay  loam,  with 
occasional  small  pebbles  present.   Texture  is 
lighter  than  B2. 
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Plot  No.  5         Bogota  Silt  Loam         Slope        1.22  percent 

Al-  0"  to  8M.   A  medium  dark  yellowish-gray  friable 
silt  loam,  with  a  weakly  granular  structure.   It 
has  a  loose  friable  consistency. 

A2-  8"  to  15".   A  light  yellowish-gray  floury  silt 
loam,  almost  ashy,  largely  single  grain  structure. 

Bl-  15"  to  18".   A  light  yellowish-gray  heavy  silt 
loam.   It  has  a  granular  structure.   With  a  brownish- 
yellow  gray  color  inside  the  granules  and  coated  with 
light  yellow  gray  on  the  outside. 

B2-  18"  to  30".   A  light  gray  clay  loam  with  splotches 
of  reddish-yellow.   It  is  rather  compact  and  plastic 
with  slow  permeability. 

C-  30"  to  40".   A  light  gray  silty  clay  loam  mottled 
with  reddish-yellow.   It  has  little  structural 
development,  and  is  lighter  in  texture  than  the  B2. 

Plot  No.  15        Elco  Silt  Loam         Slope        12.11  percent  (Eroded) 

A-  0"  to  2".   A  light  brownish-yellow  silt  loam.   A 
tendency  toward  lamination.   A  small  amount  of  clay 
is  present.   (This  layer  consists  of  a  mixture  of 
surface  soil  and  subsoil.) 

Bl-  2"  to  20".   A  light  brownish-yellow  silty  clay 
loam,  lightly  coated  and  heavily  splotched  with  gray. 
Subangular  aggregates. 

C-  20"  to  40".   A  yellowish-gray  heavy  silt  loam, 
mottled  with  brownish-yellow,  becomes  a  friable  silt 
loam  below  30".   Structureless. 

The  volume  of  storage  available  at  the  start  of  the  run  was  computed 
as  the  difference  between  total  porosity  and  soil  moisture  in  the  0"  to  21" 
depth.   Generally  the  restricting  horizon  in  these  soils  is  not  deeper  than 
21  inches.   The  volume  of  unoccupied  porosity  above  the  restricting  horizon 
at  the  start  of  the  run,  {So),   was  computed  from  antecedent  moisture 
records  and  total  porosity  determinations. 

By  inspection  of  the  last  two  columns  in  table  1,  it  is  evident  that 
the  available  porosity  is  not  always  filled  before  the  rate  of  infiltration 
becomes  constant.   The  extent  to  which  available  porosity  is  exhausted 
appears  to  be  a  function  of  the  type  of  vegetation  present.   The  volume  of 
infiltration,  which  occurred  prior  to  a  constant  rate  is  compared  with  the 
volume  of  available  porosity  for  various  crops  in  figure  2.   Plots  5  and  6 
of  figure  2  were  originally  alfalfa,  but  after  the  first  inf iltrometer  run, 
crabgrass  took  over  to  a  point  where  it  had  to  be  recognized. 


The  relationships  expressed  in  figure  2  are  not  without  precedent. 
Various  investigators  (1,  3)  have  found  that  cropping  and  depletions  in 
organic  matter  have  an  adverse  affect  on  soil  porosity.   Generally,  reduc- 
tions in  total  porosity  are  in  the  vicinity  of  10  percent  but  reductions  in 
non-capillary  porosity  are  often  more  than  50  percent.   Non-capillary 
porosity  is  the  network  of  channels  distributing  infiltrated  water  to  points 
from  whence  it  moves  by  the  slower  process  of  capillarity.   Reductions  in  the 
size,  number  or  efficiency  of  these  aqueducts  cause  rapidly  infiltrating 
water  to  bypass  a  corresponding  portion  of  the  capillary  porosity.   This 
reduces  the  portion  of  available  storage  that  will  be  filled  during  a  few 
hours  of  time  (during  the  course  of  a  rainstorm  or  irrigation  application). 

Vegetation  is  tested  here  not  for  predicting  available  porosity,  but  as 
an  index  of  the  extent  to  which  available  storage  will  be  utilized  before 
the  rate  of  infiltration  becomes  constant.   As  a  quantitative  measurement  of 
this  index,  the  basal  area,  i.e.,  the  percent  of  ground  surface  area  occupied 
by  roots  or  stems,  was  determined  for  each  plot.   Determinations  were  not 
made  for  each  run,  so  they  cannot  be  applied  individually;  however,  the 
average  density  for  each  type  of  vegetation  is  presented  in  table  2.   As 
observed  in  figure  3,  the  average  percentage  of  basal  area  is  closely 
correlated  with  the  average  percentage  of  available  porosity  filled  before 
the  rate  of  infiltration  becomes  constant. 

It  appears  that  the  potential  volume  (Fp)  of  infiltration  to  the  time 
of  constant  rate  can  be  estimated  by  applying  a  percentage  (k)  based  upon 
vegetation  to  the  initially  available  porosity  (S)  above  the  restricting 
horizon  in  the  soil.   This  procedure  provides  a  finite  value  for  reasonable 
comparison  with  known  features  of  a  soil,  and  it  places  a  ceiling  on  the 
error  of  estimate. 


CAPACITY  RATES  OF  INFILTRATION 

The  potential  volume  of  initial  infiltration  (F  @  fc)  is  of  little  or  no 
value  without  some  knowledge  of  the  rates  at  which  this  volume  will  accrue. 
Logically,  the  capacity  rate  of  infiltration  is  a  function  of  the  unoccupied 
porosity  at  a  particular  time.   Therefore,  the  status  of  remaining  potential 
infiltration  (Fp)  was  computed  in  table  1  for  comparison  with  the  rate  of 
infiltration.   Plottings  of  rate  versus  remaining  potential  volume,  on 
logarithmic  coordinates,  gave  a  curved-line  relationship,  but  subtraction  of 
the  constant  (f  )  gave  a  straight-line  relationship.   This  relationship  is 
shown  in  figure  4  for  various  types  of  vegetative  cover.   Infiltration  on 
plots  numbers  11  and  12,  Bluegrass,  was  limited  by  the  rate  of  rainfall. 
Consequently,  the  relationship  is  distorted  for  the  higher  rates  on  these 
plots  and  the  curve  was  drawn  in  with  this  in  mind.   The  general  relationship 
for  all  plots  is  expressed  by  the  formula: 

f-fc  =  aFp 
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TABLE  2. --Average  basal  areas  of  vegetation  on  inf iltrometer  plots 


Vegetation 


Number  in 
average 


Average 
basal  area 
in  percent 


Alfalfa  -  no  grass 

Crabgrass  and  alfalfa 

Weeds 

Bluegrass 

Timothy  and  lespedeza 


12 
9 
4 

11 
5 


18 
62 
10 
97 
43 
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Wherein:    f  =  capacity  rate  of  infiltration  in  inches/hour. 

fc  =  constant  rate  of  infiltration  in  inches /hour. 

F  =  potential  infiltration  to  time  of  constant  rate  and 
a_  and  n  are  constants  for  a  given  soil -vegetation 
complex. 

The  value  of  n  was  essentially  equal  to  1.387  for  all  4  curves.   The 
value  of  a  ranges  from  0.25  to  0.80  with  no  apparent  association  with 
measured  or  observed  soil  or  vegetation  influences.   Since  there  is  no 
evidence  telling  which  factors  control  the  value  of  a,  the  average  of 
a  =  0.62  was  used  together  with  the  n  =  1.387  to  represent  the  average 
relationship  illustrated  in  figure  5. 


APPLICATION  OF  RESULTS 

Pending  basis  for  further  refinement,  practical  application  can  be 
made  of  figure  5  as  follows:   (a)  Estimate  the  final  rate  of  infiltration 
from  the  soil  array  (2),  (b)  add  it  to  the  f-f_  versus  Fp  relationship  at 
points  selected  to  define  a  curved-line  relationship  of  f  to  F  ,  and  (c) 
apply  these  rates  toward  the  depletion  of  the  estimated  volume  of  the 
infiltration  potential. 

Computations  involved  in  the  inter -depletion  of  rates  and  potential 
volume  of  storage  are  illustrated  for  three  plots  in  table  3.   Although  the 
value  of  Fp  was  known  in  this  case  from  the  inf iltrometer  run,  it  is 
generally  expressed  as  kSr,  i.e.,  the  product  of  the  vegetative  factor  and 
the  available  porosity  remaining  at  any  given  time.   Actually  only  the 
initial  value  of  kSr  need  be  determined.   Subsequent  values  can  be  computed 
for  assumed  increments  of  infiltration  (  A F) .   A  corresponding  rate  of 
infiltration  (f)  is  computed  by  the  equation: 

n   *„,   !-387 
f  =  0.62kSr      +  fc 

Wherein:   f  =  rate  of  infiltration  in  inches  per  hour. 

k  =  vegetative  factor  from  figure  3. 

S  =  available  porosity  as  depleted  by  infiltrated  volumes, 
in  inches . 

fc  =  final  constant  rate  of  infiltration  in  inches  per  hour. 

The  average  rate  of  infiltration  is  applied  to  each  increment  of 
volume  to  derive  the  time  increment,  (  At).   The  curves  plotted  from 
data  computed  in  table  3  compare  quite  well  (figure  6)  with  those  observed 
with  the  inf iltrometer  on  plots  1  and  7. 
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TABLE  3. --Reproduction  of  infiltration  curves  observed  with  inf iltrometer 


kSr 

A  F 

f 

fa 

A  t 

T 

F 

(inches) 

(inches) 

(in./hr.) 

(in./hr. 

) 

(hours) 

(hours) 

(inches) 

Plot  No.  1,  Alfalfa,  run 

of  7-17- 

•40, 

kS  =  1.73,  fc  =  0.15 

o      '   c 

1.73 

0 

1.58 

0 

0 

1.50 

.23 

1.24 

1.41 

0.16       .16 

.23 

1.00 

.50 

.77 

1.00 

.50       .66 

.73 

.50 

.50 

.38 

.58 

.86      1.52 

1.23 

0 

.50 

.15 

.26 

1.92      3.44 

1.73 

Plot 

No.  7,  Lespedeza  &  timothy, 

run  of 

9-18 

-40,  kSQ  =  2.04,  fc  = 

=  0.36 

2.04 

0 

2.01 

0 

0 

1.50 

.54 

1.45 

1.73 

0.31       .31 

.54 

1.0 

.50 

.98 

1.22 

.41       .72 

1.04 

.5 

.50 

.60 

.79 

.63      1.35 

1.54 

0 

.50 

.36 

.48 
.36 

1.04      2.39 

2.04 

Plot  No.  llj 

Bluegrass,  run 

of  9-30-40, 

kSQ  =  5.72,  fc  =  0. 

L0 

5.72 

7.10 

0 

5.0 

0.72 

5.90 

6.50 

0.11       .11 

0.72 

4.0 

1.0 

4.35 

5.12 

.20       .31 

1.72 

3.0 

1.0 

2.95 

3.65 

.27       .58 

2.72 

2.0 

1.0 

1.73 

2.34 

.43      1.01 

3.72 

1.0 

1.0 

.73 

1.23 

.61      1.82 

4.72 

0 

1.0 

.10 

.42 
.10 

2.38      4.20 

5.72 

Key:    kSr  =  potential  infiltration  to  time  of  constant  rate, 
A  F  =  increment  of  infiltration. 

f  =  rate  of  infiltration  associated  with  kSr. 
fa  =  average  rate  of  infiltration. 
At  =  time  increment  =  F/f„ 
T  =  accumulated  time. 
F  *  accumulated  infiltration. 
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Figure  7  portrays  a  situation  wherein  the  computed  curve  of  infiltration 
capacity  may  be  of  material  assistance  in  analyzing  infiltrometer  data  to 
obtain  estimates  of  infiltration  capacity  during  periods  when  the  rate  of 
rainfall  is  less  than  the  infiltration  capacity  of  the  soil.   By  matching 
the  computed  rate  curve  to  those  portions  of  the  observed  curve  during  the 
later  periods  when  rainfall  was  known  to  exceed  the  infiltration  capacity, 
the  computed  curve  provides  a  projection  of  what  the  capacity  rate  would  have 
been  during  the  initial  portion  of  the  run  if  rainfall  had  been  adequate. 

The  type  F  infiltrometer,  using  a  6-  by  12-foot  plot,  has  previously 
been  demonstrated  to  give  results  comparable  to  natural  rainfall  (8).   As  a 
further  test  of  the  applicability  of  principles  derived  from  analyses  of  these 
data  to  field  scale  hydrology,  curves  of  infiltration  capacity  were  computed 
for  three  soils  for  which  adequate  data  were  available  (4,  5).   Pertinent 
characteristics  of  Kirvin  fine  sandy  loam,  Shelby  silt  loam,  and  Clinton  silt 
loam  are  given  in  table  4.   Profile  descriptions  of  these  soils  indicate  the 
restricting  horizon.   The  final  constant  rate  of  infiltration  would  logically 
be  associated  with  this  layer,'  hence,  available  porosity  was  estimated  for 
only  those  horizons  lying  above  it. 

Since  antecedent  moisture  was  not  known,  it  was  not  possible  to  check 
individual  storms.   Therefore,  the  infiltration  capacity  was  computed  for 
the  two  extreme  conditions,  saturated  and  dry.   Under  saturated  conditions 
infiltration  capacity  is  assumed  equal  to  the  final  constant  rate.   The 
constant  rates  in  table  4  were  derived  from  analyses  of  rainfall  and  runoff 
data.   However,  they  do  not  differ  greatly  from  those  estimated  for  these 
soils  from  Musgrave's  soil  array  (6).   Available  porosity  for  the  dry 
condition  was  estimated,  in  table  4,  as  the  difference  between  total  porosity 
and  the  moisture  content  at  the  wilting  point.   Porosity  below  the  wilting 
point  was  considered  to  be  available  only  through  capillary  movement, 
probably  at  a  rate  something  less  than  fc. 

Computations  of  the  infiltration  capacity  curves  for  the  dry  condition 
of  these  soils  are  given  in  table  5.   Available  porosity  derived  in  table  4 
was  modified  by  the  pertinent  vegetation  factor  k  determined  from  figure  3 
to  estimate  the  potential  volume  of  infiltration  to  the  time  of  constant  rate. 
Woods,  on  the  Kirvin  soil,  were  considered  equal  to  or  better  than  bluegrass; 
hence,  k  was  assumed  equal  to  100  percent.   Rotation  hay,  on  the  Clinton 
soil,  was  considered  comparable  to  lespedeza  and  timothy,  i.e.,  k  =  45  per- 
cent.  With  k,  SQ,  and  fc  determined,  the  computations  follow  the  same 
procedure  as  was  outlined  for  table  3.   Curves  derived  in  table  5  are  plotted 
in  figure  8. 

Records  of  rainfall  and  runoff  from  these  soils  at  Tyler,  Texas  (9), 
Bethany,  Missouri  (10),  and  LaCrosse,  Wisconsin  (2),  respectively,  were 
analyzed  to  estimate  infiltration  on  watersheds  of  several  acres  area.   The 
points  plotted  in  figure  8  represent  observed  differences  between  rainfall 
and  runoff  plotted  versus  the  total  hours  of  opportunity,  i.e.,  hours  of 
rainfall.   The  computed  extremes  essentially  bracket  the  observed  values. 
Points  lying  in  between  presumably  represent  antecedent  moisture  conditions 
intermediate  to  saturation  and  the  wilting  point. 
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TABLE  5. --Computation  of  capacity  infiltration  starting  at  wilting 
point  soil  moisture;  by  formula,  f  =  62kSr    '  +  fc 


ksr 

(inches) 


AF 
(inches) 


(in./hr.) 


(in./hr.) 


At 
(hours) 


T 
(hours) 


F 
(inches) 


Kirvin  fine  sandy  loam,  plot  GS-3,  Tyler,  Tex. ;  woods,  k  =  100  percent 


Sn  =   3.35  inches,  f  =  0.17  in/hr 


3.35 

0 

3.47 

0 

0 

0 

3.00 

.35 

3.02 

3.24 

.11 

.11 

.35 

2.50 

.50 

2.37 

2.70 

.19 

.30 

.85 

2.00- 

.50 

1.80 

2.08 

.24 

.54 

1.35 

1.50 

.50 

1.26 

1.53 

.33 

.87 

1.85 

1.00 

.50 

.79 

1.02 

.49 

1.36 

2.35 

0.50 

.50 

.40 

.60 

.84 

2.20 

2.85 

0 

.50 

.17 

.28 

1.78 

3.98 

3.35 

Shelby  silt  loam,  plot  pa-B,  Bethany,  Mo.;  Bluegrass, 
k  =  100  percent,  So  =  2.58  inches,  fc  =  0.03  in/hr 


2.58 

0 

2.40 

0 

0 

2.00 

.58 

1.66 

2.03 

0.29 

.29 

.58 

1.50 

.50 

1.12 

1.39 

.36 

.65 

1.08 

1.00 

.50 

.65 

.88 

.57 

1.22 

1.58 

.50 

.50 

.24 

.44 

1.14 

2.36 

2.08 

0 

.50 

.03 

.14 

3.57 

5.93 

2.58 

Clinton  silt  loam,  unterraced  cultivated  watershed,  LaCrosse,  Wis.; 
Rotation  Hay,,  k  est.  45%,  S0  =  6.32  inches,  f  =  0.16  in/hr 


2.84 

2.80 

0 

0 

2.50 

0.34 

2.36 

2.58 

0.14 

.14 

.34 

2.25 

.25 

2.06 

2.21 

.11 

.25 

.59 

2.00 

.25 

1.78 

1.92 

.13 

.38 

.84 

1.75 

.25 

1.50 

1.64 

.15 

.53 

1.09 

1.50 

.25 

1.25 

1.38 

.18 

.71 

1.34 

1.00 

.50 

.78 

1.02 

.49 

1.20 

1.84 

0.50 

.50 

.40 

.59 

.85 

2.05 

2.34 

0 

.50 

.16 

.28 

1.79 

4.84 

2.84 

See  key  of  table  3  for  symbols. 
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Figure  4. —  Relationships  between  capacity  rate  and  potential 
volume  of  infiltration  before  the  rate  becomes 
constant. 
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Figure  6.  —  Computed  (dashed  line)  and  observed  (solid 
line)  capacity  infiltration  on  infiltrometer 
plots  1  &  1 ,    Edwardsville,  111.   (See  foot- 
notes table  1  for  symbols.) 
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Figure  7. .-Computed  capacity  infiltration  (dashed  line)  and  observed 

infiltration  as  limited  by  rainfall  (solid  line),  inf iltrometer 
plot  No.  11,  Bluegrass,  Edwardsville,  111.   (See  footnotes  of 
table  1  for  symbols.) 
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Kirvin  f  ine  sandy  loam, woods, Tyler, Tex. 
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Shelby  silt  loam,Bluegrass, Bethany,  Mo. 
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Clinton  silt  loam,  Rotation  Hay,  LaCrosse,  Wise. 
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Hours  from  start  of  rain 

Figure  8. .-Computed  upper  and  lower  limits  (solid  lines)  of  capacity 
infiltration  compared  to  observed  amounts  (plotted  points) 
on  three  soils.   The  dashed  line  is  the  computed  upper 
limit  of  the  rate  of  infiltration. 
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